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Nonlinear  ultrasonic  harmonic  generation  to  an  input  harmonic  wave  amplitude  has 
been  simulated  first  in  ID  and  then  subsequently  in  2D  using  an  Mass  Spring  Lattice 
Model  (MSLM).  This  is  augment  to  the  earlier  developed  Finite  Difference  Time 
Domain  (FDTD)  model  in  this  project.  The  MSLM  model  is  expected  to  have  improved 
application  for  the  modeling  of  the  Non-Linear  behavior  of  the  ultrasonic  wave  propaga¬ 
tion  in  isotropic  and  anisotropic  materials  systems.  The  1-D  MSLM  model  was 
developed  and  the  NLU  behavior  has  been  simulated  and  verified  using  previously 
reported  data.  The  2-D  MSLM  model  has  also  been  developed  and  verified  for  Linear 
ultrasonic  wave  propagation  using  commercial  FEM  package.  The  incorporation  of  the 
non-linear  parameters  in  the  2-D  model  is  currently  underway. 
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1.1  INTRODUCTION 


This  report  deals  with  the  simulation  of  the  finite  amplitude  ultrasonic  wave 
propagation  in  materials  with  accumulated  no  nonlinearity.  The  wave  propagation  is 
simulated  by  solving  the  wave  equation  with  the  constitutive  behavior  of  the  nonlinear 
medium.  It  is  proved  that  an  asymmetric  stress  strain  relationship  of  the  material  result¬ 
ing  from  asymmetric  motion  of  the  dislocations  is  required  for  the  generation  of  the 
second  harmonic  and  the  static  displacement  component  while  the  third  harmonic  is 
generated  even  if  the  stress  strain  relation  is  symmetric.  The  dependencies  of  the 
generated  static  displacement  and  the  second  harmonic  on  various  input  parameters  are 
evaluated  and  are  compared  with  the  existing  literature.  The  simulation  is  carried  out 
using  a  mass  spring  lattice  model  based  finite  difference  time  domain  approach. 


1.2  MASS  SPRING  LATTICE  MODEL  -ID 

The  Mass  Spring  Lattice  model  is  an  invaluable  tool  for  simulating  wave  propa¬ 
gation  and  has  been  successfully  used  for  visualizing  waves  propagating  through 
complex  media.  The  MSLM  model  for  simulating  a  finite  amplitude  ultrasonic  wave 
propagation  through  a  material  is  based  on  the  following  assumptions.  (Holland  2002) 

(a)  The  specimen  is  a  collection  of  n  mass  points  each  of  mass  m 

(b)  The  mass  points  are  connected  by  classical  nonlinear  springs 

(c)  The  force  displacement  relation  of  the  classical  nonlinear  spring  till  the  third  or¬ 
der  is  given  by 
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(1) 


1  ,  1  3 

Fi  =klXi+-k2Xi  +~k3  Xi 


where  F-,  is  the  force  and  Xi  is  the  elongation/compression  in  the  spring. 
The  following  are  the  terms  used  in  the  derivation 


m  =  mass  of  each  mass  particle 

u[  =  displacement  of  the  particle  at  the  z-th  node  at  time  1. 

At  =  time  step 

Ax=  distance  between  two  consecutive  mass  points 
p  =  density  of  the  material 

E  =  Modulus  of  elasticity  of  the  material  (Second  order  elastic  constant) 
fi  =  Second  order  nonlinearity  parameter  of  the  material. 

A  =  area  of  cross  section  (assumed  uniform) 
y  =  Third  order  nonlinearity  parameter  of  the  material. 

Applying  Newton’s  Second  Law  of  motion  for  the  z-th  mass  we  would  have 


m  d-^±  =  (Mf+1  - u\ )  +  i k2 (u‘+1  - u\f  +  i k 3 (u‘M  - u\ )3 ]  - [kx (u\  - u'l  A ) 

at  2  3 

+  ^k2  O'  -  «Z- 1  )2  +  ^  k3  0/  -  «M  )3 1 


(2) 


i.e. 


3 


u ,+At  —2ut+u'At  1 

m  — - — 1 —  =  [k{  (ui+ 1  -  2 u,  +  mm  )  +  -  k 2  (m,+1  -  2m,  +  mm  )(m,+i  -  mm  ) 


(Ar)2 


+  -^3((m'+i  -2m'  +m'_i)((m'+1  -m')2  +(m'+1  -m')(m'  -m'_,) 
+  (m'-m'_1)2] 


(3) 


r+Af 


The  above  equation  (3)  has  only  one  term  '  referring  to  the  future  time  and  hence 
can  be  iterated  in  time  to  yield  the  displacement  matrix  u(i,t)  Re-arranging  the  above 
equation  we  have 

«r=2«;  [*,<«;*. -m -2»; 

777  2 

+  T"  ^3  ((M,+l  —  2m;  +  M;  ,  )((M,+|  —77,  )  +  (m,+1  —  M  ■  )(m,.  —  M;  j  )  +  (77,.  —  M,-  j  )  ] 

3  ••  (4) 

The  above  equation  (4)  can  be  used  to  evaluate  the  displacements  at  all  the  nodes  given 
the  initial  and  boundary  conditions  Making  clx  and  dt  tend  to  zero  in  the  (3.3)  and 
comparing  it  with  the  wave  equation  in  continuous  form  i.e. 


d2u  d2u  8  m 

^  =  £^a-2/?&-3rl 


du 

\8xj 


) 


(5) 


and  using  m=pAAx,  provides  the  following. 


k  =im 

1  Ar’  2  Ar2 


,  k , 


3  EyA 
Ar3 


The  values  of  Ax  and  At  are  chosen  to  avoid  dispersion  (as  suggested  by  Holland, 2002). 
A  Hanning  windowed  sinusoidal  pulse  was  used  as  input  to  the  simulations  conducted 
using  the  above  formulation.  A  Matlab®  code  was  utilized  to  simulate  the  wave  propa¬ 
gation  using  the  FDTD  equations  discussed  in  the  previous  section. 
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The  simulations  are  validated  against  Cantrell’s  (1984,  1987). theory  for  generation  of 


static  displacement  component 


1.2.1  Comparison  with  Cantrell’s  theory: 

Yost  and  Cantrell  (1987)  developed  an  expression  for  the  static  strain  generated  when  a 
wave  propagates  through  a  material  with  accumulated  nonlinearities 


The  equation  reads  as 


du 

dx 


1  du  f5 


C  du^ 


cn  dt  2c 


o  J 


(6) 


Where  x  is  along  the  direction  of  propagation  of  wave  u(x,t)  denotes  the  displacement  at 


du 

the  coordinate  x  at  time  1.  The  static  strain  component  is  given  by  the  time  average  of  — 


Taking  time  average  on  both  the  sides  of  equation  (6)  we  get 


1  /  du\ 
c0  \dt  / 


/(du  V 

\l  dt , 


(7) 


The  right  hand  side  of  the  above  equation  (7)  is  estimated  by  taking  the  averages  of 


from  the  simulated  u(x,t)  profiles  and  then  evaluating  the  integral  in 


the  above  equation  (7).  The  obtained  values  are  compared  with  the  static  displacement 
obtained  from  the  simulated  displacement  wave’s  Fast  Fourier  Transform.  This  was 
performed  for  different  values  of  frequency,  input  amplitude,  nonlinearity  parameter  of 
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the  material  and  propagation  distance.  One  such  comparison  is  shown  in  Table  1.  The 
results  show  excellent  agreement  between  the  published  data  and  the  simulated  data 
using  the  FDTD  MSLM  model  reported  here.  The  maximum  difference  between  the  two 
results  observed  only  in  the  second  decimal  and  can  be  attributed  to  the  numerical  errors 
during  the  computation. 

The  simulations  are  further  verified  by  comparing  variation  of  the  generated  harmonic 
with  the  distance  of  propagation,  input  amplitude  and  the  frequency  of  the  input  wave.  It 
has  been  found  to  comply  with  the  well-known  relation 

A2  gc  1 3A~k~x  /on 
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(C) 


Fig  1  shows  the  variation  of  the  second  harmonic  with  the  various  input  parameters-  the 
nonlinearity  parameter,  input  amplitude,  frequency  of  the  wave  and  the  distance  of 
propagation. 


1.3  ASYMMETRY  AND  THE  STATIC  DISPLACEMENT 
COMPONENT  GENERATION 


To  explore  the  dependence  of  the  static  displacement  component  and  the  second  har¬ 
monic  on  the  stress  strain  relation  of  the  material,  the  above  simulations  are  carried  out 
for  the  following  two  cases 

I.  /?  =  0,y  ±  0( symmetric  stress  strain  relationship) 
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II.  fi  ±  0,y=  0  (asymmetric  stress  relationship) 


Case  I: 

In  this  case  the  stress-strain  relationship  of  the  material  is  taken  in  the  following  form 

c 7  =  E(s-ys 3)  (39) 

Where  y  is  the  third  order  nonlinearity  of  the  material. 

Fig  2  shows  a  plot  of  the  simulated  time  domain  signal  of  the  wave  obtained  after  a 
distance  of  propagation  of  25  mm.  If  the  signal  is  compared  with  the  input  wave  it  is 
observed  that  the  symmetry  of  the  wave  is  restored.  In  other  words,  the  positive  half 
cycle  and  the  negative  half  cycle  distort  in  a  similar  way.  The  stress  strain  plot  for  this 
case  is  as  shown  in  Fig  3.  It  may  be  seen  that  the  stress  strain  curve  has  point  symmetry 
about  the  origin.  Thus  it  can  be  concluded  that  such  a  symmetric  stress  strain  relation¬ 
ship  of  the  damaged  material  implies  that  the  distorted  signal  restores  it  symmetry  after 
propagation. 
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Fig  2.  Wave  after  a  propagation  of  25  mm  for  Case  I. 

To  obtain  the  frequency  spectrum,  the  FFT  routines  in  Matlab  were  employed  on  the 
time  domain  signals,  for  different  distances  of  wave  propagation.  Fig  4  shows  the  FFTs 
of  the  wave  after  different  distances  of  propagation.  It  can  be  seen  that  no  static  dis¬ 
placement  component  is  generated  in  this  case,  and  only  the  odd  harmonics  are 
generated. 
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Fig  3  A  symmetric  stress  strain  curve 


Case  II: 


In  this  case  the  stress -strain  relationship  of  the  material  is  taken  in  the  following  form 


cj  =  E(£-Ps2)  (10) 

Where  fi  is  the  second  order  nonlinearity  of  the  material. 
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Fig  4  .  FFT  of  the  propagated  wave  for  Case  I  after  different  distances  of  propagation 
showing  no  generated  static  displacement  component. 

Fig  5  shows  a  plot  of  the  simulated  time  domain  signal  of  the  wave  obtained  after  a 
distance  of  propagation  of  25  mm.  If  the  transmitted  signal  is  compared  with  the  input 
signal,  an  asymmetry  induced  into  the  wave  could  be  observed  in  the  time  domain,  i.e. 
the  positive  half  cycle  and  the  negative  half  cycle  don’t  distort  the  same  way.  The 
stress-strain  plot  for  this  case  is  as  shown  in  Fig  6.  It  may  be  seen  that  the  corresponding 
stress-strain  curve  is  also  asymmetric  for  this  case. 
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Fig  7  shows  the  FFTs  of  the  wave  after  different  distances  of  propagation  for  Case  II 
material.  It  can  be  seen  that  a  static  displacement  component  increases  with  the  distance 
of  propagation.  Also  both  even  and  odd  harmonics  were  generated,  as  expected. 
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Fig  5.  Wave  after  a  propagation  of  25  mm  for  Case  II. 

Oruganti  et  al  (2007)  have  shown  that  an  asymmetry  in  the  stress-strain  curve  of  the 
material  which  arises  due  to  the  asymmetric  dislocation  motion  is  required  to  generate  a 
second  harmonic.  Similarly,  it  may  be  concluded  here  that  such  an  asymmetry  is  also 
required  to  generate  a  static  displacement  component.  It  may  be  further  noted  that  a 
symmetric  dislocation  motion  would  exist  only  when  there  are  isolated  dislocations  with 
are  initially  pinned  by  point  defects  and  are  initially  not  bowed.  Such  a  condition  would 
exist  for  example  in  an  annealed  material  which  has  very  low  dislocation  density.  For 
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materials  such  as  those  deformed  by  fatigue,  plastic  deformation  or  super  plastic  defor¬ 
mation  dislocation  motion  and  hence  an  asymmetric  dislocation  motion  would  always 
exist  and  would  contribute  to  the  generation  of  the  second  harmonic.  Further  some  of  the 
materials  processed  using  new  processes  would  have  ultra-refined  grains  and  hence  there 
would  be  an  additional  stress  caused  by  the  dislocations  in  the  cell  walls  which  form 
grain  sub-boundaries.  These  would  generate  additional  image  stresses  on  the  dislocation 
in  additional  to  the  dislocations  accumulated  within  the  cell.  These  image  stresses  would 
further  increase  the  asymmetry  in  dislocation  motion. 


Fig  6.  A  typical  asymmetric  stress-strain  curve 
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1.3.1  Dependence  of  the  static  displacement  component  on  various  input 
parameters 

The  dependence  of  the  static  displacement  component  on  the  following  input  parameters 
was  explored 

(a)  The  nonlinearity  parameter  of  the  material 

Fig  8  shows  the  variation  of  the  static  displacement  component  with  the  nonlinearity 
parameter  [i  (usually  measured  using  the  amplitude  of  the  2nd  harmonic  that  is  generated 
when  the  ultrasonic  wave  passes  through  a  non-linear  material)  of  the  material  for 
different  distances  of  propagation  for  co=5MHz  mdAi=10nm.  It  was  again  observed  that 
the  static  displacement  component  varies  linearly  with  the  nonlinearity  parameter  [>  of 
the  material  for  all  the  distances  of  propagation. 
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3.00E-06 


Fig  7.  FFT  of  the  propagated  wave  for  Case  II  after  different  distances  of  propagation 
showing  a  generated  static  displacement  component. 

(b)  The  frequency  of  the  input  wave 

Fig  9  shows  the  variation  of  the  static  displacement  component  with  the  frequency  of  the 
input  wave  for  different  distances  of  propagation  for  f=16  and  Ai=10nm. 

From  the  least  square  fit  it  can  be  seen  that  the  static  displacement  component  varies  as 
the  square  of  the  frequency  of  the  input  wave  for  all  the  distances  of  propagation. 
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Fig  8.  Variation  of  the  static  displacement  component  with  the  nonlinearity  paramter  for 
different  distances  of  propagation  for  co-  5 MHz,A]  =  10 run . 


(c)  The  amplitude  of  the  input  wave 

Fig  10  shows  the  variation  of  the  static  displacement  component  with  the  amplitude  of 
the  input  wave  for  different  distances  of  propagation  for  co=5MHz  and  f=16  .  From  least 
square  fit  it  can  be  seen  that  the  static  displacement  component  varies  as  the  square  of 
the  amplitude  of  the  input  wave  for  all  the  distances  of  propagation. 


(d)  The  distance  of  propagation 

Fig  11  shows  the  variation  of  the  static  displacement  component  with  the  distance  of 
propagation  for  various  input  amplitudes  for  co=5MHz  and  [i= 1 6.  It  was  observed  that 
the  static  displacement  component  varies  linearly  with  the  distance  of  propagation. 
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Fig  9.  Variation  of  the  static  displacement  component  with  frequency  for  different 
distances  of  propagation  for  /?  =  16,  A1  =  10 run . 
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Fig  10.  Variation  of  the  static  displacement  component  with  the  amplitude  of  the  input 
wave  for  different  distances  of  propagation  for  an  input  wave  of  frequency 

co  =  5MHz,j3  =  16. 
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Fig  11.  Variation  of  the  static  displacement  component  with  the  distance  of  propagation 
for  different  input  amplitudes  for  an  input  wave  of  frequency  co  =  5 MHz,  ft  =  16 


1.4  2D  ULTRASONIC  WAVE  PROPAGATION  MODEL  FOR 
MSLM 


Mass  Spring  lattice  Model  is  one  among  the  numerical  simulation  technique  for  model¬ 
ing,  simulating  and  visualizing  elastic  wave  phenomena  by  discretizing  the  material  into 
a  collection  of  masses  interconnected  by  springs.  In  the  model  two  types  of  springs, 
linear  and  torsional  are  present.  This  is  well  illustrated  in  Fig.  12. 
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Fig  12.  Schematic  of  a  2D  MSLM  Model. 


1.5  THE  2D  MSLM  MODEL 

The  equation  of  motion  for  the  centre  of  mass  in  the  x  and  y  direction  can  be  written  as 
+^ij  -2«5)/(Ar)=  Fx+kx{u-+lj  +u-_lj-2u.j)/h2  +k2(u-+lJ+1  +u-+lj_1  +<1J+1 


+  uti  i- 1  -4wM/2/i2  +k2/3(ii 


k  .  k  .  k  .  k 
i+\,j+\  +  Mi+l,y-l  +  +  Ui-\,j-\ 


-4u*  j)/2h2  +k2(v 


2  \  y  i+lj+l 


~vlu-i  ~ vi-i,j+i  +  -4u.j)/2h2  +k2/3(-v.+l  j+1  +v*u_1  +v*u+1  -v,*^  -4u.  ])l2h7 


-2v.j)/h2  +k2(v.+lJ+l  +vf+1J_j  +vf_1J+1 


P(v“  +v,j‘  -2v‘y)/(A(2)=  F,  +i3(v‘u  +v‘1J  -2v.  j)/h2  +v‘lj_, 

+  vf-ij-i  -4vfj)/2h~  +k2fi(vkM  j+l  +v.+lJ_l+v._lJ+1  +v._lJ_l  -4v)j)t2h2  +k2(u.+lJ+1  ~u.+lJ_ 


+  m*1m  -4m*  A/2h2 +k2fi(-u 


j+l.j+l  +  Mi+1,7-1  +  M  (-1,7+1  Ui-l,j 


<i.m-4</)/2^2 


where  wf -and  vf  ■  denote  the  x  and  y  displacements  on  the  mass  particle  at  the  position 
(i,j)  and  at  the  time  t  =  kAt , where  At  is  the  time  step \kl,  k2,  /cl, and  a  are  the  stiffnesses. 
/?  can  be  expressed  as  ft  =  a/(2h2 )  . 
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The  wave  equations  in  plain  strain  transverse  isotropic  elastic  medium  is 
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0  u 
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0  u 
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d  u 
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dr 
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ay- 


+  (C,2  +  C33) 


d~u 

ydxdy  j 


+  C 


f  rs2  \ 

O  V 


>■ 


33 


dx1 


+  F, 


(12) 


where  u  and  v  are  the  displacement  in  the  x  and  y  direction  respectively, 
p  is  the  mass  density  of  the  medium,  CL  are  the  elements  of  stiffness  matrix.  For  iso¬ 
tropic  materials  CL  are  expressed  in  terms  of  the  Lame’s  constant  A  and  p  as 
C[  j  =  C2 2  =  A  -\-  p  C33  —  p  snd  C|  2  —  A . 


After  discretizing  the  above  equation  using  center  difference  method  and  as  the  LHS  of 
two  equations  are  equal,  RHS  can  be  equated,  then  we  get 


sr  +r  c  -c 

=  C, ,  — C33  ,k3  = — — - —  ,  k3-C22-C33,  and  /?  = — — - —  for  transversely 

4  i  C\2 

isotropic  material  and  kl=k3  =  A  +  p  ,k2  =  ^  +  and  ft  =  — — —  for  isotropic  media 

4  A  +  3p 

From  the  first  set  of  equations  the  displacement  for  the  (k+1)  th  time  step  can  be  deter¬ 
mined  from  the  displacement  values  of  k  thand  (  k-1)  thtime  step  as  follows. 


(«u  =~uij  +2ui.j  +F,/p(At2)+kl(u-+lJ  +u._hj  -2m; kj)/p(At2)h2  +k2(u-+lj+1  +u-+lj_l  +u-_h 
+  uUj- 1  ~4ui, j)/ p(&2)lh2  +  k2fi(u-+lj+1  +u.+lJ_1+u._hj+1+u._l  j_1  -4ul)l p(\t2)lh2  +k2(y[ 
~vLj- 1  ~vi-u+i  +vmj-i  ~4u. pi p(\t2)lh2  +k2P(-v-+l  j+l  +vf+1J_i  +  vi-ij+i  -vll  j_1-4u-  j)/p(At 


3+1 


+1,3+1 
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(v“  =-v*j'  +2v5  =FJ/p(Af!)+M(vLJ  +vLj  -2v5)/p(A(J>J  +  ^(v,tw  +v,‘;1J_1  +v*. 
+  Vwj-,  -4v‘,)/p(A(5)2/i!  +vLj_,  +v‘,.w  +v,V.  -4v'j)/p(At2)2A2  +*,(«*,. 

-“Lj.i  -4«,‘j)//’(a/2)2As  +k2P(-ullajtt  +u‘tlJ_,  +u‘_tJt, -4h‘j)/p(a/’)2/i 


7+1 

_  * 

7+1  Mi+1.7'-1 

2 


where  h  is  the  space  step  (grid  size  )  and  df  is  the  ime  step. 


Element  size  and  time  step  determination 


To  accurately  model  the  wave,  there  should  be  at  least  8  elements  within  the  minimum 
wavelength.  Minimum  wavelength  is  for  shear  wave  and  hence  for  the  present  case  the 
grid  size  is  taken  as 


h  =  -fr  and  2lrtn  =  T-  , 

16  / 


where  A  is  the  wave  length  Cs  is  the  shear  velocity  in  solid  and  /  is  the  frequency.  The 
shear  wave  velocity  can  be  obtained  from  the  material  parameters  E  and  Poisson’s  ratio 
as 


s  2p(l  +  v)' 

For  a  stable  wave  to  propagate  through  the  explicit  numerical  scheme  the  Courant 

number,  ?  where  C  is  the  longitudinal  wave  velocity.  Hence  the  time  increment 

h 

is  taken  as  dt  - 


h 

42  c  ’ 


where  C  = 


E{  1-v) 

Pi  l  +  v)(l-2v) 
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As  per  Hyunjune  and  Younghoon  (2000),  after  a  critical  value  of  1.6937xl04  for  — 

h 

the  numerical  scheme  is  not  stable.  For  our  particular  case  it  is  1.3964x10  4.  So  stability 
condition  is  met  with  the  present  MSLM. 


Fig.  13.  The  selection  of  critical  values  using  as  per  Hyunjune  and  Younghoon  (2000) 


1.6  2D  MSLM  SIMULATION  RESULTS  USING  MATLAB 

Material  selected  for  the  simulation  is  steel  with  density  7800  kg/m“  ,  Poisson’s  ratio  0.3 
and  Young’s  modulus  200  MPa.  The  material  is  considered  as  a  collection  of  masses  and 
spring  and  a  MATLAB  code  was  written  for  finding  the  displacement  in  the  x  and  y 
direction  for  the  equation  of  motion. 

Figure  .14a,  shows  the  wave  propagation  when  a  single  cycle  sine  pulse  was  given  at  the 
centre  and  imposing  rigid  boundary  condition.  Both  longitudinal  and  shear  waves  can  be 
seen.  When  the  excitation  was  given  at  the  surface  Fig.  14b,  surface  waves  are  absent 
with  rigid  boundary  condition. 


22 


4.3853e-006 


Fig.l4a.Source  at  the  centre  Fig.  14b  Source  at  the  surface 

For  the  proper  visualization  of  free  surface  condition,  equations  are  written  separately 
for  the  nodes  along  the  boundary  and  at  the  corners.  Along  the  boundary  masses  and 
springs  are  halved  and  at  the  comers  only  !4  th  mass  is  taken.  After  incorporating  the 
above  condition  surface  wave  could  also  be  seen  in  addition  to  longitudinal  and  shear 
waves  as  shown  in  Fig  15. 

5.8471  e-006 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

50  100  150  200 


Fig.  15.  The  visualization  of  both  L  and  S  waves  using  MSLM  2D  model. 


Later  code  was  modified  by,  giving  frequency,  number  of  cycles  and  dimensions  of  the 
specimen  as  input.  The  results  for  a  3  cycle  100  kHz  hanning  windowed  pulse  is  shown 
below  for  different  time  intervals  after  the  initial  pulse  was  applied  on  the  surface.. 
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Fig.  16.  The  2D  visualization  of  wave  propagation  in  a  2D  isotropic  media  at  different 
time  intervals  using  the  MSLM  model. 


‘A’  scans  at  different  points  were  also  plotted  as  shown  in  Fig.  17  at  discrete  points  in 
the  2D  media. 
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Fig. 17  A  scans  at  different  points. (a). At  the  source  point(l,128)  (b)  At  (50,128)  (c)At 
(80,128)  and  (d)At  (128,128) 

1.7  COMPARISON  WITH  COMMERCIAL  FEM  SOFTWARE 
(ABAQUS). 

Results  obtained  with  MSLM  are  compared  with  the  results  obtained  from  Abaqus  in 
which  the  same  source  was  given  as  the  input.  The  wave  propagation  simiulation  plots 
are  provide  in  the  form  of  color  images  at  different  time  intervals  in  Fig.  18  The  plots 
with  GREY  background  are  resutls  obtained  using  ABAQUS  while  the  plots  with 
WHITE  background  represent  the  MSLM  model  results  in  Fig.  19. 
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(a)  ABAQUS  FEM  MODEL 


(b)  MSLM  MODEL 


Fig.  18.  The  simulation  of  the  wave  propagation  using  ABAQUS  FEM  model  for  a 
points  source  in  a  2D  domain  at  2  different  time  intervals,  using  (a)  ABAQUS  and  (b) 
MSLM  Model  developed  in  MATLAB. 
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(a)  At  the  source  point 


(b)  At  the  right  corner 


(c)At  the  right  edge  middle  point 

Fig.  19.  The  A-scan  signals  at  different  points  as  obtained  from  ABAQUS  and  MSLM 
models. 
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It  is  concluded  from  these  results  that  the  2D  MSLM  model  developed  using 
MATLAB,  developed  in  this  project  compares  well  with  the  2D  ABAQUS  FEM  model 
that  is  commercially  available.  Hence,  the  validity  of  the  MSLM  model  for  linear 
ultrasonic  wave  propagation  is  validated  herewith.  The  advantage  of  the  MSLM  model, 
over  commercially  available  FEM  models  includes  (a)  The  ability  to  model  the  mechan¬ 
ics  of  the  wave  propagation  at  micron  scales,  and  (b)  The  ability  to  introduce  non-linear 
behavior  of  the  springs,  which  can  be  used  to  simulate  NLU  parameters  for  the  current 
studies. 


1.8  FUTURE  WORK  IN  2D  MSLM  MODELLING 


It  is  planned  to  introduce  non-linear  behavior  of  the  individual  components  of  the 
MSLM  model  in  order  to  simulate  NLU  behavior,  as  was  earlier  demonstrated  in  the  ID 
case.  Once  this  is  developed,  the  3D  version  will  be  developed  for  the  improved  under¬ 
standing  of  the  NLU  behavior  in  materials. 
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